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Abstract 

Eleven morphometric and one meristic character in 15 parthenogenetic Artemia populations have been studied by 
using discriminant and cluster analysis as well as scanning electron microscopy. Discriminant analysis revealed 
five main groups of morphological patterns: (i) the coastal Chinese populations together with a population from 
Kazakhstan, (ii) the inland Chinese salt lake populations, (iii) the Greek populations, (iv) one African population 
from Namibia and (v) a Chinese population from Xuyu (Jiangsu province). Cluster analysis was not always in 
agreement with discriminant analysis and these results are discussed. 


Introduction 

The genus Artemia comprises bisexual species, which 
are found on all continents except Antarctica, and 
parthenogenetic populations, which are endemic in 
Europe, Asia and Australia. Artemia is known for 
two phenomena which are rare in the animal king¬ 
dom: obligatory parthenogenesis and polyploidy. Poly¬ 
ploidy is common among parthenogenetic populations 
which are either diploid (2 n = 42), triploid (3 n = 63), 
tetraploid (4 n = 84) or pentaploid (5n = 105) with 
diploidy and tetraploidy being the most frequently 
observed (for a review see Barigozzi, 1974; Abreu- 
Grobois, 1987). Mixed ploidy levels often occur in nat¬ 
ural populations (Abreu-Grobois & Beardmore, 1982; 
1991; Abatzopoulos et al., 1986; Pilla, 1992; Zhang 
& King, 1993). Amat (1980) and Hontoria & Amat 
(1992a) studying parthenogenetic populations from the 


Western Mediterranean basin demonstrated that ploidy 
level affects the morphology of Artemia. 

Ionic composition of the habitat can produce eco¬ 
logical isolation (Bowen et al., 1985, 1988) and can 
result in morphological differences (Hontoria & Amat, 
1992b). 

In this paper we study parthenogenetic populations 
living in geographically isolated habitats and we try 
to analyze their morphological characteristics using 
multivariate methods. 


Materials and methods 

Fifteen parthenogenetic populations were evaluated in 
this study. Table 1 summarizes the origin of each pop¬ 
ulation and the abbreviations that are used. For the 
Chinese populations we refer also to the code number 
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Table 1. List of the studied populations, Artemia Reference Center (ARC) cyst bank code 
number and abbreviations. For the Chinese populations the code numbers proposed by Xin 
et al. (1994) are given as well 


Population 

ARC cyst 
bank number 

Code number 

of Chinese 
populations 

Abbreviation 

Megalon Embolon, Greece 

1279 


MEM 

Citros, Greece 

1280 


CIT 

Polychnitos, Greece 

1281 


POL 

Kalloni, Greece 

1282 


KAL 

Aibi Lake, Xinjiang, P. R. China 

1198 

11010891 

AIL 

Balikun Lake, Xinjiang, P. R. China 

1235 

11030991 

BAL 

Kazakhstan (unknown source)* 

1060 


CAT 

Chengkou, Shandong, P. R. China 

1210 

04010991 

CHK 

Dongjiagou, Liaoning, P. R. China 

1216 

01030991 

DOJ 

Hangu, Tianjin, P. R. China 

1212 

03010991 

HAN 

Swakopmund, Namibia 

1186 


NAM 

Pulandian, Liaoning, P. R. China 

1217 

01040991 

PUL 

Xuyu, Jiangsu, P. R. China 

1283 

05010992 

xuw 

Huanghua, Hebei, P R. China 

1233 

02041091 

HUA 

Daban Lake, Xinjiang, P. R. China 

1197 

11020891 

DAL 


* Sample received from Catvis bv, Netherlands. 



Figure I. Scatterplot of the first three canonical discriminant functions (group centroids) resulting from the discriminant analysis when using 
the origin of each population as separation factor. 


proposed by Xin et al. (1994) in order to better describe 
the different strains. 

The methodology for studying the morphology of 
the adults, the culture conditions, the methodology for 
scanning electron microscopy as well as the statistical 
analysis have been described in detail in a previous 


study (Triantaphyllidis et al., 1997). For each popula¬ 
tion at least 28 individuals have been measured. Sur¬ 
vival ranged between 56 ± 12% and 86 ± 9% for the 
populations studied. 








Table 2. Mean values (standard deviation in parenthesis) of morphometric and meristic characters of different parthenogenetic Artemia populations, n = number of animals analysed. 
Abbreviations of populations in Table 1. A: total length, B: abdominal length, C: length from the third abdominal segment to the end of the abdomen, D: length of the eighth abdominal 
segment, E: width of third abdominal segment, F: length of furca, G: width of head, H: length of first antenna, I: distance between eyes, J: diameter of complex eye, K: number of setae 
on the left branch of the furca, L: number of setae on the right branch of the furca, M: width of ovisac. Wilks’ lambda (A) is given by the equation 1 — rj 2 , where r} 1 is the ratio of 
the between-groups sum of squares to the total sum of squares and represents the proportion of the total variance attributable to differences among the groups (Norusis, 1993) 
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Results 

Chi-square and Kolmogorov-Smirnov tests showed 
that the transformed morphological characters fit to the 
normal distribution (P > 0.05). Two different analy¬ 
sis of variance tests were performed for each of the 
thirteen characters listed in Table 2. The first was a 
standard one-way analysis of variance (Anova), where 
variances are assumed to be homogeneous (Sokal & 
Rohlf, 1981). However, Levene’s test (Norusis, 1993) 
showed that the variances were significantly differ¬ 
ent and since logarithmic or square root transforma¬ 
tions did not eliminate the problem, we conducted an 
approximate test of equality of means according to the 
procedure of Games & Howell (1976). 

The mean values of morphometric and meristic 
characters are presented in Table 2. Anova resulted 
in F ratios, the significance of which revealed high¬ 
ly statistical significant differences (P < 0.00005). 
The variables that present the larger F ratios are the 
length from the third abdominal segment to the end of 
the abdomen and the abdominal length. Table 3 sum¬ 
marizes the significant differences of the means after 
applying the test of Games & Howell (1976). 

Discriminant analysis based on the origin of each 
population as a separation criterion resulted in 13 
canonical discriminant functions after 26 steps. The 
first nine discriminant functions are highly statistically 
significant (P < 0.0001) and obtain a cumulative sep¬ 
aration percentage of 99.42% while the last four are 
not statistically significant (P > 0.05). The first three 
functions that appear in Table 4 resulted in a separation 
percentage of 78.23%. In Table 4 the unstandardized 
and standardized canonical discriminant function coef¬ 
ficients are also presented. Using morphometric and 
meristic characters each population can be classified 
correctly in one of the fifteen groups with an overall 
accuracy of 93.20% (Table 5). 

Figure 1 depicts the plot of the discriminant analy¬ 
sis, while Figure 2 summarizes the results of hierar¬ 
chical cluster analysis using the single linkage method 
(nearest neighbor technique) in a dendrogram. NAM 
population is clustered separately from the other popu¬ 
lations suggesting that it has unique morphological fea¬ 
tures. The Greek populations CIT, KAL and POL are 
grouped very close while MEM appears to be clustered 
a little bit further from the ‘Greek core’. XUW popu¬ 
lation is discriminated from all the other populations. 
The Chinese populations CHK, HAN, PUL and DOJ 
(all located in coastal areas) exhibited great morpho¬ 
logical similarities. Close to them appeared the CAT 



Figure 2. Dendrogram of hierarchical cluster analysis using the 
single linkage method to join clusters (nearest neighbor technique). 


population which inhabits an inland salt lake. BAL, 
DAL and ATI, are close in the three dimensional space 
of Figure 1, suggesting that these inland salt lake popu¬ 
lations from Xinjiang province (P. R. China) exhibited 
morphological similarities. AIL is closer to DAL while 
BAL is further discriminated by the first two functions. 
Although in Figure 1 discriminant analysis separates 
the Chinese inland-lake populations BAL, DAL and 
ATT, and groups them closely, cluster analysis classi¬ 
fies BAL population separately, close to the coastal 
habitat group (see Figure 2). 

Summarizing the above data, discriminant analysis 
revealed five main groups of morphological patterns. 

(i) The coastal Chinese populations CHK, PUL, HAN, 
DOJ and HUA together with the CAT population 
from Kazakhstan, 

(ii) the inland Chinese salt lake populations BAL, DAL 
and AIL, 

(iii) the Greek polyploid populations CIT, KAL, POL 
and MEM, 

(iv) the African NAM population and 

(v) the Chinese XUW population. 

Scanning electron microscopy of the head and 
uterus showed the same difficulties of quantification as 
in the bisexual animals (Triantaphyllidis et al., 1995a). 
The scanning electron micrographs of Figure 3 depict 
an example of the differences observed in the morphol¬ 
ogy of the furca. 


Discussion 

The coastal Chinese populations are all diploid (Pilla, 
1992; Triantaphyllidis et al., unpublished data) with 
the exception of HAN for which few tetraploid indi¬ 
viduals appeared. Because the habitat of these popu¬ 
lations is coastal they all experience an environment 
with similar ionic composition (chloride is the prevail¬ 
ing anion). CAT population, from Kazakhstan, seem 
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Table 3. Significant differences for the mean values of morphometric and meristic characters of parthenogenetic females as determined by 
the approximate test of equality of means of Games & Howell (1976). Populations that share the same letter (a. b, ..., h) per row, are not 
significantly different (a=0.05). The codes (A, B,.... M) for the various characters are the same as in Table 2. The abbreviations of populations 
can be found in Table 1 



MEM 

CIT 

POL 

KAL 

AIL 

BAL 

CAT 

CHK 

DOJ 

HAN 

NAM 

PUL 

XUW HT JA 

DAL 

A 

e,f,g,h 

e.fig,h 

g.h 

b,c,d,e 

e.f.g 

b.c.d.e 

c,d,e,f 

b 

b,c,d 

b,c 

a 

b,c,d,e 

e,f,g,l 

i c,d,e,f 

h 

B 

e,f,g 

e,f,g,h 

e,f,g 

b,c,d,e 

c,d,e,f,g 

c,d,e,f,g b,c,d,e 

b 

b,c 

b 

a 

b,c,d 

e,fig 

b,c,d,e,: 

f h 

C 

d,e,f 

e,f,g 

e,f 

b,c,d,e 

e,f 

d,e,f 

b,c,d,e 

b,c 

b,c 

b 

a 

b,c,d 

e,f 

c,d,e,f 

g 

D 

b,c,d,e,f,g d,e,f,g,l 

i fig.h 

b,c,d,e,f,g d,e,f,g,h 

d,e,fig.h b,c,d 

b,c 

b,c,d,e 

b 

a 

b,c,d,e,1 

f d,e,fig d,e,f,g,h g.h 

E 

f,g,h 

b,c,d,e,i 

f d,e,f,g,l 

i f,g,h 

c,d,e,f 

a,b 

c,d,e,f,g,h 

a 

b,c,d,e,l 

F a,b,c,d 

c,d,e,f,g a,b,c 

fig.h 

b,c,d,e 

b,c,d,e,f 

F 

b,c,d 

b,c,d 

b,c,d 

a,b,c 

a,b 

a 

e 

b,c,d 

d 

d 

e 

b,c,d 

b,c,d 

b.c.d 

b,c,d 

G 

h 

e,f,g 

f,g,h 

b,c,d,e,f 

e,f,g,h 

a,b,c,d 

b,c,d,e 

a,b,c,d,e 

a,b,c 

a,b 

a 

a,b,c 

e.fig 

d.e.fig 

b,c,d,e,f 

H 

e,f 

d,e 

e,f 

a,b,c,d 

b,c,d 

a,b 

a,b,c 

a,b,c,d 

b,c,d 

a,b 

a 

b,c,d 

f 

e 

b,c,d 

I 

f,g 

c,d,e 

d,e,f 

b,c,d 

c,d,e 

b 

b,c 

b,c,d 

b 

a,b 

a 

b 

g 

c,d,e,f 

b,c,d,e 

J 

% 

c,d,e,f 

c,d,e 

a,b,c,d 

d,e,f,g 

a,b,c 

c,d,e,f,g 

c,d,e,f,g 

c,d 

a,b,c,d 

a 

a,b,c 

h 

d,e.fig 

a,b,c,d 

K 

a,b 

a 

a,b 

a,b,c 

a,b,c 

c,d,e 

f 

c,d,e,f 

e,f 

c,d,e,f 

f 

c,d,e,f 

a,b,c 

b,c,d,e 

b,c,d 

L 

a 

a 

a 

a,b 

a,b,c 

b,c,d,e 

f 

c,d,e,f 

e,f 

d,e,f 

f 

e,f 

a,b,c 

c,d,e,f 

b,c,d,e 

M 

g.h 

b,c,d,e 

d,e,f,g 

d,e,f,g,h 

b,c,d 

a,b 

b,c,d,e 

a 

b,c,d,e 

b,c,d 

d,e,f 

b,c,d 

h 

b,c 

b,c,d,e 


Table 4. Results of the discriminant analysis. The classification was based on the origin of each population 


Variable 

Unstandardized canonical discriminant 

function coefficients* 

Standardized canonical discriminant 

function coefficients* 

1 

2 

3 

1 

2 

3 

Total length 

20.5503828 

-1.9919531 

-0.9182448 

0.59319 

-0.05750 

-0.02651 

Abdominal length 

-8.7883584 

2.8049734 

3.3432993 

-0.31134 

0.09937 

0.11844 

Length up to third ab.s. 

26.6589978 

-4.1719837 

2.1435532 

0.96804 

-0.15149 

0.07784 

Length of the 8th abd.s 

1.9909037 

-2.0349200 

-3.1557132 

0.09085 

-0.09286 

-0.14400 

Width of 3rd abd. segm. 

-8.4306974 

-1.5479552 

-8.6719410 

-0.27365 

-0.05025 

' -0.28148 

Length of furca 

-7.3195161 

2.2719251 

-0.1865445 

-0.48763 

0.15136 

-0.01243 

Width of head 

-4.2482949 

-16.0326828 

-29.3804486 

-0.11865 

-0.44778 

-0.82057 

Length of 1st antenna 

-5.6265499 

7.7990300 

0.1409685 

-0.18562 

0.25729 

0.00465 

Distance between eyes 

7.6838676 

26.1122522 

24.1504402 

0.22133 

0.75214 

0.69563 

Diameter of eyes 

-3.4419944 

6.4880032 

18.3382265 

-0.11871 

0.22376 

0.63246 

No. of setae left branch 

0.5978158 

-0.5033924 

1.3742926 

0.18628 

-0.15686 

0.42822 

No. of setae right branch 

0.1405275 

-0.7587268 

1.2522186 

0.04300 

-0.23218 

0.38320 

Width of the ovisac 

-11.9591987 

7.5841354 

-2.5473703 

-0.54556 

0.34598 

-0.11621 

Constant 

-40.8168018 

0.2035977 

-7.2071804 






Eigenvalue** 

Percentage of 
variance 

Cumulative 

percentage 

Canonical 

correlation 

Wilks’ 

lambda*** 

Chi-square 

DF 

P 

Function 1 

5.8768 

53.78 

53.78 

0.9244 

0.003982 

2436.967 

182 

< 0.0001 

Function 2 

1.4922 

13.66 

67.44 

0.7738 

0.027383 

1586.651 

156 

<0.0001 

Function 3 

1.1786 

10.79 

78.23 

0.7355 

0.068242 

1183.948 

132 

<0.0001 


* The unstandardized coefficients are the multipliers of the variables when they are expressed in the original units, while the standardized 
coefficients are used when the variables are standardized to a mean of 0 and a standard deviation of 1, just like in multiple regression 
(Norusis, 1993). 


** Eigenvalue is the ratio of the between groups to within groups sums of squares. Large eigenvalues are associated with ‘good’ 
functions (Norusis, 1993). 

*** Wilks’ lambda (A) is given by the equation A = 1 - tj 2 , where ij 1 is the ratio of the between-groups sum of squares to the total 
sum of squares and represents the proportion of the total variance attributable to differences among the groups (Norusis, 1993). 
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Figure 3. Scanning electron micrographs of the differences observed in the length of the furca and number of setae. Plate A: NAM population, 
plate B: CIT population. 
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Table 5. Classification results of discriminant analysis showing the percentage of individuals classified in each group. The diagonal 
elements are the number of cases classified correctly into the groups and serve as an indicator of the effectiveness of the discriminant 
analysis. The percent of ‘grouped’ cases correctly classified is 93.20%. The abbreviations of populations can be found in Table 1 


Actual No of Predicted Group Membership (%) 


group 

Cases 

MEM 

CIT 

POL 

KAL 

AIL 

BAL 

CAT 

CHK 

DOJ 

HAN 

NAM 

PUL 

XUW 

HUA 

DAL 

MEM 

31 

93.5 

0.0 

3.2 

0.0 

3.2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CIT 

31 

0.0 

90.3 

6.5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3.2 

POL 

31 

3.2 

6.5 

83.9 

3.2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3.2 

KAL 

33 

9.1 

0.0 

3.0 

87.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

AIL 

30 

3.3 

0.0 

0.0 

0.0 

90.0 

3.3 

0.0 

0.0 

0.0 

0.0 

0.0 

3.3 

0.0 

0.0 

0.0 

BAL 

30 

0.0 

0.0 

0.0 

0.0 

0.0 

100.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CAT 

30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CHK 

30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

93.3 

0.0 

0.0 

0.0 

3.3 

0.0 

3.3 

0.0 

DOJ 

30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

93.3 

0.0 

0.0 

6.7 

0.0 

0.0 

0.0 

HAN 

30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3.3 

93.3 

0.0 

3.3 

0.0 

0.0 

0.0 

NAM 

32 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100.0 

0.0 

0.0 

0.0 

0.0 

PUL 

30 

0.0 

0.0 

0.0 

3.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

93.3 

0.0 

0.0 

3.3 

XUW 

30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100.0 

0.0 

0.0 

HUA 

30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3.3 

0.0 

0.0 

10.0 

0.0 

83.3 

0.0 

DAL 

28 

0.0 

0.0 

0.0 

0.0 

0.0 

3.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

96.4 


to exhibit great morphological similarities to the pre¬ 
vious populations. It is known that CAT population is 
diploid (Pilla, 1992) but we lack information of the ion¬ 
ic composition and of the exact geographical location 
of this population. Hontoria & Amat (1992a) showed 
that parthenogenetic diploid and tetraploid popula¬ 
tions, from the Western Mediterranean basin, can be 
thoroughly recognized by their morphological charac¬ 
teristics. The same ploidy level of the populations that 
form the first group seemed to result in similar mor¬ 
phological patterns. 

BAL, DAL and AIL populations are in geographi¬ 
cally close inland salt lakes where sulphate is the pre¬ 
vailing anion. They presented different ploidy levels: 
AIL is diploid, BAL tetraploid and DAL a mixture 
of diploids and tetraploids (Pilla, 1992; K. Thomas, 
personal communication). The differences in the ploi- 
dies could be the reason why DAL population appears 
between AIL and BAL. 

The third group is comprised by the Greek popula¬ 
tions. They are polyploid (Abatzopoulos et al., 1986; 
Triantaphyllidis et al., 1993) and they all exist in 
coastal saltworks. Although MEM population appears 
to diverge from the ‘core’ of the Greek cluster, the 
Greek populations are clearly discriminated from all 
the other parthenogenetic populations. 

NAM population is a mixture of different ploidies 
(Barigozzi, 1986). However, the studied batch is main¬ 
ly diploid (Triantaphyllidis et al., unpublished data) 


and exhibits some characteristics that allow 100% dis¬ 
crimination from all the other populations studied here. 
The telson is very long (the longest among the popu¬ 
lations studied here) and the abdomen is very short, 
characteristics that allow the NAM population to form 
a group alone, far discriminated from all the other 
groups. 

Due to lack of information about the ploidy level of 
the XUW population it is rather difficult to understand 
the reasons that contributed to such distinct discrimi¬ 
nation. Discriminant analysis classifies it with 100% 
prediction and this suggests that the overall morphol¬ 
ogy of this population differs not only from the other 
Chinese populations but from all the other populations 
considered in this study. Cluster analysis groups it close 
to the Greek populations and this could possibly mean 
that it is a polyploid population. Further study needs to 
be done on this population. 

Geographically-isolated populations, even if they 
share the same ploidy level, do not necessarily present 
similar morphological patterns. For instance, although 
BAL population and the four Greek populations are 
all tetraploid, they do not cluster together; this may 
be due to the fact that BAL population inhabits a sul¬ 
phate environment while the Greek populations live in 
coastal salinas where chloride is the prevailing anion. 
Despite their similarities in their ploidy level the lat¬ 
ter populations experience environments with different 
ionic composition and it is possible that natural selec- 
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tion favours genotypes (clones) which are better adapt¬ 
ed to that particular environment. This selection could 
eventually lead to slight morphological changes that 
hecome apparent after applying multivariate statistical 
methods. It is known that ecological isolation exists 
in Artemia and that there are differences in tolerance 
of anion concentrations between populations (Bowen 
et al., 1985, 1988). Ecological isolation or environ¬ 
ments with different ionic composition can lead to mor¬ 
phological differentiation (Hontoria & Amat, 1992b). 
Thus, ploidy level is possible to dictate a morpholog¬ 
ical pattern only in geographically-close populations 
provided that the environmental or physical traits of 
their ecosystems are also quite identical. As more pop¬ 
ulations are made available for study, we see that the 
possible morphological repertoires are enriched. 

Parthenogenetic populations can be discriminated 
from females of bisexual populations as has already 
been demonstrated by Triantaphyllidis et al. (1995) 
for an A, franciscana and a parthenogenetic popula¬ 
tion from Tanggu (Tianjin, People’s Republic of Chi¬ 
na). By applying discriminant analysis the morpholog¬ 
ical differences between parthenogenetic and bisexual 
females allowed allocation to the correct species as 
of day 14 with almost 100% accuracy and in salin¬ 
ities ranging from 35 to 180 ppt (Triantaphyllidis et 
al., 1995b). The same resulted from merging the data 
of this study with the data of the earlier study of Tri¬ 
antaphyllidis et al. (1995). In all cases parthenogenet¬ 
ic females could be discriminated from the bisexual 
females with 100% accuracy. 

The results from the discriminant analysis present¬ 
ed in Figure 1 seem to contradict the results of cluster 
analysis presented in Figure 2, i.e. BAL population 
in discriminant analysis seems to cluster closer to the 
coastal Chinese group rather than to the Chinese inland 
lake group. This partial conflict between the results of 
discriminant and cluster analysis was not observed in 
the case of bisexual populations (Triantaphyllidis et 
al., 1997). The potential problems of cluster analysis 
occur at every level of the strategy (Pimentel, 1979). A 
given set of cluster analyses can be subjected to more 
than one interpretation. The hierarchical feature limits 
the method to extracting a one dimensional portrayal 
of relationships among individuals or populations in 
our case. Farris (1977) has demonstrated that the sin¬ 
gle dimension extracted by cluster analysis need not 
be the most parsimonious representation of the oper¬ 
ational taxonomic units (OTU’s). Whenever relation¬ 
ships among OTU’s require more dimensions for clear 
protrayal, cluster analysis will be imprecise. Since 


different strategies can result in drastically different 
dendrograms, various interpretations of a given study 
are possible. Choosing the ‘best’ strategy should be 
approached very carefully and must be accomplished 
on a non-statistical basis (Pimentel, 1979); different 
researchers, most probably, would select different den¬ 
drograms as being most biologically meaningful. ‘In 
spite of the above difficulties, there is little doubt that 
cluster analysis can and has produced useful results. 
This would imply that in given circumstances (as in 
the case of bisexual populations described in a previ¬ 
ous study) the potential pitfalls either do not exist, or 
methodology is sufficiently robust to circumvent the 
pitfalls’ (Pimentel, 1979). 


Conclusions 

Application of multivariate methods in the morpholog¬ 
ic study of Artemia populations seems to be a useful 
tool. If information of the chromosome number, ionic 
composition of the water as well as climatic conditions 
will be combined with genetic data (i.e. isozymes) 
then we can obtain a very precise picture of the stud¬ 
ied populations. In agreement with Hontoria & Amat 
(1992a,b) great care should be given to the standard¬ 
ization of the culture conditions and we repeat here the 
proposal for using standard laboratory conditions for 
all the studies that will take place in the future. 
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Abstract A scanning electron microscopy (SEM) study 
of bisexual Artemia populations revealed that popula¬ 
tions representing the species A. franciscana, A. per- 
similis, A. urmiana , A. sinica and a recently described 
species from Kazakhstan have a pair of spine-like out¬ 
growths at the basal parts of their penes, whereas pop¬ 
ulations from southern Europe and North Africa (i.e. 
Mediterranean populations) lack these spine-like out¬ 
growths. Allozyme and DNA polymorphisms, detected 
by allozyme starch gel electrophoresis and AFLP fin¬ 
gerprinting, respectively, suggested conspecificity of the 
studied populations from the broader Mediterranean 
basin. Male specimens from the collection of the Natural 
History Museum of London (UK) of the extinct A. sa- 
lina population from Lymington lack spine-like out¬ 
growths at the basal parts of the penes. This finding, 
based on a taxonomic character which is quite reliable, 
suggests conspecificity of A. salina from Lymington and 
the present bisexual Artemia populations from the 
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Mediterranean basin, grouped under the binomen 
A. tunisiana. 


Introduction 

The brine shrimp Artemia (Crustacea: Branchiopoda) is 
a well-studied organism, however, taxonomists are still 
puzzled about the evolution and the phylogenetic rela¬ 
tionship of the populations that comprise the genus (for 
a review of the confusing names in the genus Artemia see 
Belk and Brtek 1995). Schlosser made the first descrip¬ 
tion of the brine shrimp in 1755 on material collected 
from solar saltworks near Lymington, England, UK 
(Kuenen and Baas-Becking 1938). Linnaeus described 
the brine shrimp as Cancer salinus in 1758 and Leach 
renamed it as A. salina in 1819 (Artom 1931). In the 
following years many populations have been identified, 
and nowadays the genus Artemia is comprised of a 
complex of bisexual species and superspecies as well as 
parthenogenetic populations with various degrees of 
ploidy (Browne and Bowen 1991). 

Artom (1905, 1906, 1907), Stella (1933), Stefani 
(1963), Stefani and Cadeddu (1967), Halfer Cervini et al. 
(1968), Piccinelli et al. (1968) and Barigozzi (1974), 
studying the Artemia populations in Italy, used the 
binomen A. salina for the bisexual populations that had 
42 chromosomes. Ever since 1910 (Barigozzi 1974) and 
even recently (Sorgeloos and Beardmore 1995), most 
scientific papers referred to all brine shrimps as A. salina 
although in the meantime it was known that two distinct 
modes of reproduction occurred (parthenogenetic or 
zygogenetic) and that new bisexual species had been 
described: A. franciscana (Kellogg 1906) in the New 
World and A. urmiana (Gunther 1890) in Lake Urmia, 
Iran. Piccinelli and Prosdocimi (1968) described a new 
species, A. persimilis, from the saltworks of San Barto¬ 
lomeo, Cagliari (Sardinia) and Hidalgo (Argentina). 
This species has 44 chromosomes and a different adult 
morphology. It was considered to be sympatric with 
A. salina. Clark and Bowen (1976) showed that the 








